Double Ionization by Strong Elliptically Polarized Laser Pulses 
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We join the tribute to Professor N.B. Delone in this memorial issue by presenting the results of 
new calculations on the effects of ellipticity on double ionization by short and strong near-optical 
laser pulses. 

PACS numbers: 32.80.Rm, 32.60.+i 



Observations of double ionization by femtosecond-scale 
near-optical laser pulses with intensities above about 
lO^^W/cm^ revealed a prominent surprising feature. The 
doubly ionized ion count, as a function of incident laser 
intensity, was observed to be significantly above the value 
predicted by the famous ADK ionization formula [ij, 
which is used to describe non-correlated tunneling es- 
cape of two electrons in sequence. The presumption is 
then natural that a strongly correlated ionization chan- 
nel has become open. The ADK formula and its prede- 
cessors 0, 111 are well matched by experimental facts for 
single ionization so it was startling to find them violated 
for high-field double ionization. Recent overviews of this 
domain of laser-atom physics are given by Agostini and 
DiMauro 3 and Becker and Rottke 

A three-step mechanism Q has been generally ac- 
cepted as responsible for highly correlated double ion- 
ization: first, one electron is freed with zero kinetic en- 
ergy by tunneling through the coulomb barrier, which 
has been tipped by the laser field but is nearly static on 
the Bohr orbit time scale (Keldysh approximation [2]); 
second, the freed electron is first pulled away and then 
driven back when the oscillating laser field reverses; and 
third, it acquires enough acceleration from the laser field 
to eject another electron in a near-core collision process. 
The same recollision process plays a key role in high har- 
monic generation Q- 

In this short note we present the first data from our 
theoretical analyses that include ellipticity of the laser 
field, suggesting this traditional view is too rigid, and 
may need modification. While the calculations them- 
selves are not of the type that Professor Delone would 
have undertaken, as we will describe, we believe he would 
have been intrigued by the results emerging. The conse- 
quences of elliptical polarization, as an influence in high- 
field time-dependent atomic phenomena, are not at all 
trivial, but early reports mentioning tunnel ionization 
for elliptical polarization [l], 0] indicate that for cases not 
too close to circular one has a simple ratio of linear to 
elliptical tunnel rates that is independent of field inten- 



sity: 



Wellip. 
Wlin. 



1 



(1) 



where e is the ellipticity. 

In Fig. [l]we show three double ionization probability 
curves as a function of e. They indicate predictions of 
double ionization for e — Q (linear polarization), e = 0.5, 
and £ = 1 (circular polarization). As one expects, the 
linear polarization curve remains a substantial fraction 
higher than the other two. However, it is intriguing that 
at the lower-intensity end of the curves all three deviate 
noticeably from the well-known smoothly and strongly 
decreasing curve typically provided by the ADK formula. 
We believe that this indicates an interesting domain for 
experimental work because such a deviation from the 
ADK curve is usually regarded as more or less impos- 
sible for finite ellipticity because the return collision in 
the three-step ionization mechanism is difficult to justify. 

Our approach, by which the curves in Fig. [T]were cal- 
culated, is heretical in the sense that it is fully classical, 
without even allowing tunneling to initiate the first elec- 
tron's escape. However, a fully classical trajectory analy- 
sis fSf] is well known for NSDI. In this alternative analysis 
two electrons can share energy inside the classical core 
potential and occasionally liberate one of them without 
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FIG. 1: Predictions for double ionization as a function of laser 
intensity for ellipticity e = 0, 0.5, and 1.0. 
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FIG. 2: A trapezoidal pulse with phase = 0. The curve of higher 
amplitude is the x-component of the field and the lower curve is 
the y-component. The ellipticity is here taken to be £ = 0.5. In 
calculations the phase in Eqn. JSj is set randomly to mimic the 
phase variation of repeated laser pulses. 

tunneling. This picture has been satisfactory in explain- 
ing at all experimental double ionization (NSDI) features 
observed under linear polarization and has been ex- 
tended to triple ionization \ld^ in good agreement with 
existing data, e.g., the relative shapes and strengths of 
momentum distributions in double ionization compared 
to triple ionization of argon [llj . 

Here we have used a two-dimensional ensemble 0, 
where the total energy of the system can be written as: 

z z 



V(a;i-a;2)2 + (yi-y2)2 + fe2' 

The first term is the total kinetic energy of the system, 
the second and third terms are soft-core Coulomb at- 
tractions Q between the nucleus and the two electrons 
and the fourth term is the soft-core Coulomb repulsion 
between the two electrons. The size of the ensemble is 
chosen between 500K and lOM, depending on different 
ellipticities, which is large enough to get a statistically 
unbiased subensemble of double ionization events. We 
have set the total energy of each 2e member of the en- 
semble at -1.3 a.u., which is close enough to the ground 
state energies of both Kr and Xe to allow fairly close 
experimental relevance. 

A many-pilot-atom method p^ . in the absence of a 
laser field, is used to generate a microcanonical zero-field 
ground state ensemble to start the interaction at t=0. A 
laser pulse as shown in Fig. [2] is used, where the linear 
two-cycle ramp has some theoretical advantages [l3| . Ap- 
proximately, it gives an iso-intensity interaction between 
field and electrons, but it sacrifices smoothing associated 
with more realistic pulses shapes. 



Once the microcanonical ensemble is generated, the 
laser field 

E{t) = Eof{t)[exSm{Ljt + (f>) + eyecos{ujt -t- 0)], (3) 

is turned on. The circular frequency lu is set to be 0.0584 
a.u., corresponding to 780 nm wavelength, consistent 
with most experiments using Ti:sapphire laser systems. 
A 10-cycle pulse with this frequency gives a pulse length 
of about 26 fs. 

The classical potential's parameter values a — 1.77 and 
6 = 0.1 serve to specify a model for the atomic core. The 
amplitude of the laser field is set to give a peak (not 
average) laser intensity of 0.6 PW. At the end of the laser 
pulse, the positions and momenta of the two electrons are 
recorded. The end-of-pulse ion counts are shown in Fig. 
[T]and the momenta of the doubly ionized ions, which are 
calculated as the sum of the momenta of the two ionized 
electrons, are shown in Fig. [3jLeft) as a distribution in 
the Px — Py plane. 

A scan of Fig. [3](Left) over the distribution of ion 
momenta in the a;-direction is a single-peak curve similar 
to what is seen with linear polarization. However, the 
y direction is more interesting because it contains the 
consequences of non-zero e. There are four rather sharp 
peaks of momenta, as shown in Fig. [31[Right). These 
are features not available previously, and they appear to 
deserve careful future attention. 

Experience has shown that experiments in the PW 
regime of intensities and multi-femtosecond pulse dura- 
tions provide detailed information about the double ion- 
ization process, particularly by examining the recoil ion 
momentum distribution. For elliptic or circular polar- 
ization, because of the dramatic drop of the probability 
of the first-freed electron to return directly to the ionic 
core, high-correlation ionization has not been detected in 
most rare gas atoms or molecules. Of course, one cannot 
automatically conclude that there is no NSDI for elliptic 
or circular polarization, and in fact characteristic knee 
signatures have been observed in the double ionization 
of Mg [l3| and several molecules even with circular po- 
larization [l5|. This is some support for the suggestion 
made above that the three curves in Fig. [T]show depar- 
tures from ADK-type predictions. 

These results are preliminary and not conclusive, but 
they represent an end-of-pulse examination of double 
ionization under ellipticities. They indicate that the 
full range of ellipticities from linear to circular polariza- 
tion can be analysed using our purely classical ensemble 
method. The novel four-band Py structure identified for e 
= 0.5 also appears in preliminary results for other values 
of e. We plan to extend these calculations, and expect 
on the basis of previous success with the classical model 
that prominent features appearing, such as the four-peak 
character of Fig. [3l will be able to be explained by a 
straightforward trajectory analysis. 
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FIG. 3: Effects of elliptical polarization for e = 0.5. In the px — Py plot (Left) one sees a structured end-of-pulse momentum distribution 
for doubly ionized ions. A scan produces the py distribution shown in the Right. 
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